How to suppress undesired synchronization 
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It is delightful to observe the emergency of synchronization in the blinking of fireflies to attract 
partners and preys. Other charming examples of synchronization can also be found in a wide 
range of phenomena such as, e.g., neurons firing, lasers cascades, chemical reactions, and opinion 
formation. However, in many situations the formation of a coherent state is not pleasant and 
should be mitigated. For example, the onset of synchronization can be the root of epileptic seizures, 
traffic congestion in communication networks, and the collapse of constructions. Here we propose 
the use of contrarians to suppress undesired synchronization. We perform a comparative study 
of different strategies, either requiring local or total knowledge of the system, and show that the 
most efficient one solely requires local information. Our results also reveal that, even when the 
distribution of neighboring interactions is narrow, significant improvement in mitigation is observed 
when contrarians sit at the highly connected elements. The same qualitative results are obtained 
for artificially generated networks as well as two real ones, namely, the Routers of the Internet and 
a neuronal network. 



I. INTRODUCTION 

In the year 2000, Londoners were presented with the 
Millennium Bridge, a futuristic footbridge that became 
the center of attention on the inauguration day. The el- 
bowing of the crowd, eager to be the first to cross it, 
forced the synchronization of walkers causing a lateral 
swing of the structure pQ. Once on this wobbly struc- 
ture, how could one avoid such uncomfortable situation? 
Abnormal synchronization is also the origin of neurolog- 
ical diseases such as epilepsy and Parkinson [2]. Brain 
pacemakers have been developed and implanted in the 
patient to discharge an electrical signal into the brain 
tissue and restore the normal activity [3HH] . But imagine 
a, still to develop, device able to interact with individ- 
ual neurons. What would be the best strategy to break 
the synchronization? A third source of inspiration can 
be found in the Internet, where several interconnected 
Routers receive and redistribute the information pack- 
ages in the network. When multiple Routers synchronize 
their delivering events, the network collapses, a dysfunc- 
tion known as TCP global synchronization. To avoid it, 
several algorithms have been developed and implemented 
in some Routers |F . What is the fraction of such proac- 
tive Routers required to avoid global synchronization? In 
social context, avoiding synchronization might represent 
a political tool to fight a charismatic leader. Consider a 
speech that inflames a crowd. Initially every individual 
claps at his/her own rhythm but rapidly a coherent clap 
emerges 0. If a set of political adversaries (contrarians) 
try to destroy the harmony, what would be the best strat- 
egy, the proper amount of contrarians, and their spatial 
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distribution in the hall? 

The Kuramoto model has extensively been used as the 
paradigm to study synchronization [SHllj. In a first at- 
tempt to address the questions raised above, we gener- 
alize this well-established model to include contrarians 
which try to suppress the emergence of global synchro- 
nization. We present a systematic study of how the syn- 
chronizability depends on the fraction of contrarian os- 
cillators for two different strategies and analyze the influ- 
ence of the topology in the mitigation process. To illus- 
trate our results, contrarian oscillators have been studied 
in silico for two real networks, namely, the Routers that 
compose the Internet |12j and the network of neurons 
of the organism C. elegans [T3l [14] . Our results suggest 
that local contrarians can be used as a powerful way to 
control synchronization, avoiding the necessity of moni- 
toring the global state. Moreover, spreading contrarians 
as hubs is also much more effective. 



II. MODEL 

The described examples are characterized by a set of 
N oscillators (walkers, neurons, Routers, or spectators), 
mutually interacting. Hereafter we take the example of 
the walkers but the model can be straightforwardly ex- 
tended to all other cases. The stepping of each walker 
i is characterized by the phase Oi(t) and its natural fre- 
quency Ui, corresponding to the stepping frequency when 
isolated. When the crowd moves, all walkers initially step 
at their natural frequency but herding (under strong cou- 
pling) rapidly leads to coherent walking [T]. In the Ku- 
ramoto model, the motion of each oscillator is described 
by a phasor e 1 ^ 4 ', where 9i(t) is the phase, and the cou- 
pling between walkers is such that the dynamics of each 
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is governed by, 

JV 

0i =uj 1 + \J2 A ij sin (Oj - °i) > (!) 

where the sum is over all other walkers (i ^ j), A is the 
coupling strength, is the walker's natural frequency, 
and A is the connectivity matrix such that Aij = 1 if 
walker i is influenced by walker j or zero otherwise. The 
collective walking can be characterized by the complex 
order parameter defined as, 

K*y* (t) = ^fy^- (t) , (2) 

3 

where the sum is over all walkers, the amplitude 
< — 1 measures the global coherence, and 

is the average phase. 

To account for contrarians we introduce a second pop- 
ulation of N c walkers also coupled with the others but fol- 
lowing a different dynamics. A contrarian k is also char- 
acterized by the phase Ok{t) and its natural frequency 
LOk- We consider two types of coupling: a mean-field 
(Model A) and a pairwise (Model B). In the mean-field 
coupling the dynamics of contrarians is governed by, 

d k =u k + \sm(*-O k -S), (3) 

where ^(t) is the average phase and 6 is a phase shift. 
In the pairwise coupling the dynamics of contrarians is 
governed by, 

N T 

9 k = uj k + A ^ A kj sin (6j - 9 k - 6) , (4) 

j=i 

where the sum is over all the Nt walkers (Nt = N + N c ). 
Hereafter, we take 5 = ir in both cases. Such phase shift 
between two walkers (k and j) would correspond to a 
walking such that when k steps with the left foot j steps 
with the right. It is noteworthy that these two mod- 
els yield different types of frustration. While in Model A 
frustration between regular and contrarian walkers is me- 
diated by the average phase, in Model B the frustration 
results from a pairwise interaction between regular and 
contrarian walkers where the former attempts to mutu- 
ally synchronize while the latter tries to dephase. 

Two models accounting for frustration in the mean- 
field Kuramoto model have recently been discussed in the 
context of a mixture of positive and negative couplings. 
Zannette [15] considered a pairwise coupling where the 
strength and sign of the interaction between two oscil- 
lators is symmetric. This model, in the limit uj k = 0, 
is equivalent to a magnetic XY model with a distribu- 
tion of couplings. Hong and Strogatz proposed a differ- 
ent scheme, where regular walkers are also solely coupled 
with the average phase and the spatial distribution 
of regular walkers is not considered [THl HZ] . This model 



is similar to the mean-field limit of Model A discussed 
here. In contrast to the model discussed here, synchro- 
nization cannot be suppressed in any of these previous 
models. 



III. RESULTS 

In the absence of contrarians, the classical Kuramoto 
model is characterized by the emergence of synchroniza- 
tion at a critical coupling A c which depends on the dis- 
tribution of natural frequencies (ui) and on the degree. 
While under weak coupling (A < A c ) the motion is in- 
coherent (r — 0), above the critical coupling a coherent 
motion emerges (r > 0). In the limit of very strong cou- 
pling (A 3> A c ) all oscillators participate in the coherent 
motion. 

The presence of contrarians can affect the coherent mo- 
tion. In Figure [I] different fractions p = N c /N of con- 
trarian oscillators are considered in the mean-field (A) 
and pairwise (B) models. While mean-field contrarians 
are not able to reduce the value of r, a fraction as small 
as 5% of pairwise contrarians is enough to significantly 
reduce the synchronizability. Further investigation shows 
that, although contrarians enable the system to desyn- 
chronize, mean-field contrarians drive the system to a po- 
larized state, where oscillators are concentrated around 
two phases: — 7r and tt. It is possible to understand this 
splitting through the analysis of the stable point (given 
by Ok = 0) for contrarians in the mean-field model, yield- 
ing 

Asin($ - e k - tt) = -w fe . (5) 

Assuming that uj k is symmetrically distributed around 
zero, this equation shows that the difference between the 
phase of contrarians and the average phase must be equal 
to tt. Thus, contrarians have a tendency towards the ex- 
tremes of the possible phases, dragging their conformist 
neighbors in the process. Hence, mean-field contrarians 
introduce differences in the dynamic behavior of oscilla- 
tors by polarizing them in two distinct phases (see Sup- 
plementary Fig. SI). 

For the pairwise (B) contrarians, the emergence of a 
coherent state is suppressed above a certain fraction of 
contrarians (Figure [2]). We notice that synchronization is 
suppressed for values of p > 0.15, as shown by the peak in 
the standard-deviation of r in the inset of Figure [2] The 
peak increases with the network size. For small values 
of p synchronization is maintained, r > 0, as conformists 
synchronize their phases with each other and the small 
fraction of contrarians dephase from their neighbors with- 
out destroying global synchronization. In this situation, 
the average phase of contrarians and conformists create a 
periodic alternating wave over time, an interesting mech- 
anism that resembles, for instance, the oscillation of pop- 
ulations of predators and preys which characterizes the 
classical Lotka-Volterra model [T5] ■ Figure [3] shows an 
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FIG. 1. Comparison between mean-field and pairwise 

coupling. Dependence of the order parameter r on the cou- 
pling strength A, for the mean-field (A) and the pairwise (B) 
couplings. Different curves stand for different fractions of 
contrarians p randomly distributed in a random graph with 
average degree equal to four. 



example where contrarians (the central layer of the net- 
works in the upper part) are in opposition to their first 
neighbors conformists, which in turn try to synchronize 
with them. The sequence of networks in the figure are 
snapshots of the oscillators and their phase over time. 
Before all conformists could match their phases with con- 
trarians, the latter already have an opposite phase. The 
periodic wave that conformists and contrarians create is 
clear in the lower part of Figure [3] which shows the aver- 
age phase of different types of oscillators. 

As the fraction of contrarians overcomes a certain 
threshold, the effect of contrarians spreads over the en- 
tire network completely suppressing global synchroniza- 
tion. This suppression is a consequence of an increasing 
fraction of contrarian/contrarian interactions, which nat- 
urally tend to be dephased, reinforcing their impact. For 




FIG. 2. Impact of pairwise contrarians on the synchro- 
nization. Order parameter r dependence on the fraction of 
contrarians p showing suppression of synchronization. Differ- 
ent curves stand for different network sizes. The inset contains 
a plot of the standard deviation of r among samples showing 
a transition around p — 0.15. A coupling strength of A = 2.0 
has been used. 



p > 0.15, neither synchronization (r = 0) nor a periodic 
wave is observed (see Supplementary Fig. S2). 



A. Contrarians as hubs 

The phenomenon of synchronization is known to result 
from the interplay between the network topology and the 
dynamics of oscillators [19 , 20 . In this section, we discuss 
the improvement in the desynchronization efficiency by 
distributing the contrarians among the oscillators (nodes) 
with higher degree and compare this strategy with the 
random distribution case discussed above. 

We start considering the case of a random graph 
(Erdos-Renyi (ER) network), characterized by a Poisson 
distribution of degree. As shown in the Supplementary 
Fig. S3, in spite of the narrow degree distribution, the 
fraction of contrarians necessary to reduce synchroniza- 
tion is reduced to one third (p = 0.05) when contrarians 
sit at the nodes with higher degree. In this case, the 
disturbing effect of contrarians occurs even for smaller 
fractions p and synchronization is efficiently suppressed. 
In scale-free networks, where degree distribution is scale 
free, and highly connected nodes are more frequent, the 
assignment of hubs as contrarians favors desynchroniza- 
tion (see Supplementary Fig. S4). 

In real networks, the presence of communities and 
other features, such as assortativity and clustering, also 
play a role in synchronization [2U [22] . Here we consider 
two real networks and show that the same behavior holds 
(see Figure |4| . The first one is the network of Routers 
in the Internet. This network is believed to have grown 
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FIG. 3. Time dependence of the phase. Upper part: Snapshot of a network of 200 oscillators at 4 different time steps 
(vertical dashed lines in the lower part). A total of 20 pairwise contrarians are displayed in the central layer. Each concentric 
layer i, from inside to outside, contains the ith neighbors of the contrarians. The color of each node represents its phase. Lower 
part: Time dependence of the average phase of contrarians (red squares), conformists (blue triangles), and the whole set of 
oscillators (black circles) showing a periodic oscillation over time. 



through the mechanism of preferential attachment, being 
characterized by a scale-free degree distribution of expo- 
nent 7 = 3.00 [23]. Moreover, it has been shown to have 
a hub dominant structure, where many hubs share low 
degree neighbors [53]. In this case, similarly to scale-free 
networks, hubs play a major role in the synchronization, 
and only 10% of them are necessary to suppress syn- 
chronization. Highly connected contrarians have a disor- 
dering effect on a great number of conformists and a few 
hubs control the phase of the entire set of oscillators. The 
second one is the neuronal network of the organism C. 
elegans, a Small- World network and the largest network 
of neurons that has been totally mapped [53]. There, 
random or degree-based distribution of contrarians sup- 
press synchronization, although the network seems to be 
more resilient to such control than ER networks. The 
presence of functional communities of neurons that are 
highly connected within themselves might be the cause 
of this resistance [26]. It is interesting to note that this 
biological system, evolved under evolutionary pressure, 
has converged to a resilient structure regarding synchro- 
nization. 



IV. DISCUSSION 

The use of agents to control the dynamics of a network 
has recently been a subject of research interest [2"T1 128j . 



In this work, we have shown that global synchroniza- 
tion can be suppressed with local agents (contrarians) 
which systematically dephase from their nearest neigh- 
bors. We show that solely local information is required 
to efficiently avoid a coherent oscillatory state. If in- 
stead, global information is considered, the set of oscilla- 
tors splits into two oscillatory states with different phases 
and a global coherent state is still possible. Analyzing 
the impact of the network topology in desynchronization 
we concluded that, when contrarians sit at the nodes 
of higher degree, the process is more efficient than in 
the random case. Even with random graphs, character- 
ized by a narrow degree distribution, the degree strategy 
reduces to one third the amount of required contrari- 
ans to suppress synchronization. We also show that the 
synchronization of real networks that present underlying 
features such as communities and dominant hubs is also 
suppressed with the use of contrarians. 

The present work is the first attempt to understand 
the interplay between the desynchronization dynamics 
and the topology when considered to mimic real systems. 
The Kuramoto model provides a standard framework to 
study synchronization, however it entails several approx- 
imations when discussing real networks. For example, so- 
cial systems are composed of adaptive agents that might 
change their strategy over time and avoid being trapped 
in a locked state. Developments build up on top of this 
work should account for further details on the coupling 
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FIG. 4. Pairwise contrarians on real networks. Fraction of contrarian oscillators randomly assigned and based on the 
degree for different networks: (a) the Routers of the Internet, (b) neurons on C. elegans. The insets are snapshots of the 
referred networks where each concentric layer i, from inside to outside, contains the ith neighbors of the contrarians. The color 
of each node represents its type: contrarians (green) and conformists (yellow). A coupling strength of A = 4.0 has been used. 



scheme and on the contrarian dynamics. For instance, 
here we have focused on the distribution of contrarians, 
but recent works have shown that rewiring interventions, 
such as swapping, adding, or removing edges, have a cru- 
cial role in the collective dynamics [29, .30;. Nevertheless, 
general conclusions can be drawn shedding light on the 
problems discussed in the introduction. For instance, as 
referred, some Routers on the Internet have a special al- 
gorithm implemented to avoid synchronization. We have 
shown that placing contrarian Routers as hubs on the 
network could optimize the fraction of proactive Routers 
necessary to prevent global synchronization. Also, mod- 
ifying the coupling mechanism between two Routers and 
the contrarian strategy implemented, it becomes possible 
to extend our work to define the best location of contrar- 
ians on the Internet. 

Regarding the development of brain pacemakers, our 
study suggests that a set of small size devices spread 
throughout the brain, and solely tracking the phase of 
neighboring neurons, would be more effective to prevent 
a seizure than monitoring the global state. After the 
mapping of a neural network [31 , and the characteriza- 
tion of the coupling dynamics between neurons, our work 
also gives helpful hints about the minimum amount and 
the optimal spatial distribution, of these active devices. 

Interesting applications could also be found in so- 
cial dynamics. Whenever a "social synchronization" is 
achieved, such as clapping after a speech, a small amount 
of influential agents can be trained to prevent this syn- 
chronization. Political opponents could be spread in the 
crowd to avoid a proper salutation simply by "dephasing" 
their claps with their close neighbors. The same method 
might be used to prevent a synchronized walk on a bridge 
where instructed actors could walk dephased from others. 
Evidently crowd behavior control is a very hard task |32j , 
but here we show that it could theoretically be achieved. 



V. METHODS 

Equations [l] |3j and [4] have been numerically solved 
using a fourth order Runge-Kutta method with discrete 
time steps St = 0.001. The complex order parameter was 
computed in the interval t € [90, 100] using the average 
value of Equation [2j For all considered cases, natural 
frequencies of oscillators have been uniformly distributed 
between —0.5 and 0.5 and initial phases have also been 
uniformly distributed between —n and n. A coupling 
strength of A = 2.0 has been used. For the Internet and 
the C. elegans, natural frequencies and initial phases have 
been distributed uniformly between —0.1 and 0.1, and 
between — f and § , respectively, and a coupling strength 
of A = 4.0 has been used. 

Network Nodes Average Degree Max Degree 
The Internet 40028 2^36 259 
C. elegans 283 17.39 115 

TABLE I. Number of nodes, average degree, and maximum 
degree for network of Routers in the Internet and the neural 
network of the C. elegans. 

Networks of oscillators have been constructed as undi- 
rected ER networks of average degree four, unless oth- 
erwise stated. The network of Internet Routers has 
been analyzed through data retrieved from the Opte 
Project that represents all the communication among 
40028 Routers on January 15th of 2005 [U]. Each node 
of this network is a Router with an associated IP ad- 
dress and the links (edges) are established between two 
IP address which have communicated at least once. The 
network of neurons was constructed through data ob- 
tained on the Worm Web website jUj and is mostly based 
on the work by Chen et al 13J. In this network, links 
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have been established whenever an interaction between 
neurons has been registered, regardless of their type or 
direction. Other measures regarding these networks are 
available in Table HI 

All results have been averaged over several samples. 
The error bars were omitted in all figures, being smaller 
than the symbols. For the Internet and the C. el- 
egans only the natural frequencies and initial phases 
change among samples. Plots in Figure [I] are con- 
structed using the average value over 100 networks of size 
N = 1000. Figure [2] shows averages of 5000 networks of 
size N = 1000, 1000 networks of size N = 5000, and 600 
networks of size N = 10000. In Figure [3j we represent 



a single network of N — 200 with p = 0.1 and A = 1.0. 
For Figure |4j plot (a) is an average over 200 initial dis- 
tributions of phases and frequencies, and plot (b) is an 
average over 10. 
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FIG. SI. Distribution of oscillator phases. Histogram of 
the phase of oscillators for t = 100. Curves represent oscilla- 
tors according to the mean- field (A) and pairwise (B) models. 
While mean-field contrarians are frozen in — tv and n, pairwise 
contrarians are uniformly distributed through all phases. Re- 
sults are averages over 1000 ER networks of 1000 oscillators 
each, where p — 0.1 and A = 1.0. 
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FIG. S2. Average phase of oscillators with different fractions of pairwise p contrarians. Time dependence of the 
average phase of oscillators for different fractions of randomly assigned contrarians. The amplitude of the wave goes to zero 
with the fraction contrarians oscillators. 




FIG. S3. Impact of pairwise contrarians assigned to 
the nodes of highest degrees on the synchronization 
of ER networks of average degree four. Order parame- 
ter r dependence on the fraction of contrarians p for different 
network sizes showing a suppression of synchronization af- 
ter the introduction of pairwise contrarians. The inset is the 
standard-deviation of r showing a transition around p = 0.05, 
much smaller than ramdomly assigned contrarians. 
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FIG. S4. Scale-free networks with contrarians. Fraction of pairwise contrarian oscillators assigned randomly and based 
on their degree to scale-free networks of different degree exponent 7, namely, 1.75, 2.5, and 3.25. 



